Formation of molecules from a Cs Bose-Einstein condensate 
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Conversion of an expanding Bose-Einstein condensate of Cs atoms to a molecular one with an 
efficiency of more than 30% was observed recently in experiments by M. Mark et al, Europhys. 
Lett. 69, 706 (2005). The theory presented here describes the experimental results. Values of 
resonance strength of 8 mG and rate coefficients for atom-molecule deactivation of 1 x 10~ n cm 3 /s 
and molecule-molecule one of 1.5 x 10~ 9 cm 3 /s are estimated by a fit of the theoretical results to the 
experimental data. Near the resonance, where the highest conversion efficiency was observed, the 
results demonstrate strong sensitivity to the magnetic field ripple and inhomogeneity. A conversion 
efficiency of about 60% is predicted by non-mean-field calculations for the densities and sweep rates 
lower than the ones used in the experiments. 

PACS numbers: 03.75.Mn, 03.75.Nt, 82.20.Xr 
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Introduction 

A molecular Bose-Einstein condensate (BEC) has been 
recently formed in experiments on atomic BEC Q, 
El [J and on quantum-degenerate Fermi gases 
The molecules have been formed by sweeping the Zee- 
man shift through a Feshbach resonance (see Ref. 0) m 
a backward direction, so that the molecular state crossed 
the atomic ones downwards. This led to the transfer of 
population from the lowest atomic state in the case of 
a BEC, or from an energy band in the case of a Fermi 
as, to the molecular state, as had been proposed in Ref. 
Assuming all the atomic population is initially in 
the BEC state, the backward sweep would have been 
ideally suitable for forming molecules, were it not for 
two loss mechanisms. The resonant molecule is gener- 
ally populated in an excited rovibrational and electronic 
state, and therefore can be deactivated by exoergic in- 
elastic collisions with atoms and other molecules (see 
Refs. @ E3, El E2) (A stabilization of molecules by 
means of coherent control has been considered in 
Ref. \14\.) In addition, during the backward sweep, some 
higher-lying non-condensate atomic states can be popu- 
lated temporarily due to molecular dissociation. These 
two effects restrict the efficiency of conversion from the 
atomic BEC to the molecular one. 

Formation of molecular BEC from degenerate Fermi 
gases, realized in experiments 0, is more efficient due 
to Pauli blocking of inelastic collisions Some pecu- 
liarities of this process have been discussed in Ref. [T^ |. 
In the case of Bose atoms the deactivation losses can be 
minimized by a reduction of the condensate density (see 
Ref. [13 )■ Such a reduction can be realized in experi- 
ments by using an expanding BEC, released from a trap 
by turning it off (see Refs. 0, H H 3). A conversion 
efficiency of more than 30%, comparable to the one ob- 
tained for Fermi gases, has been achieved by this way in 
Cs experiments |2(. 

The present work provides a theoretical description of 
molecular formation in an expanding BEC with applica- 
tions to the Cs experiments 2] , including an explanation 



of the high conversion efficiency. Some of the results 
obtained below can be satisfactory derived by a mean- 
field theory. However, in general, a description of cer- 
tain processes, such as the spontaneous dissociation of 
the molecular BEC into non-condensate entangled atom 
pairs discussed in Sec. Q] below requires the use of a non- 
mean-field theory. Several theoretical methods are avail- 
able for this purpose. One such method is based on a 
numerical solution of stochastic differential equations in 
the positive- P representation, as used in the present con- 
text in Ref. Another method is the Hartree-Fock- 
Bogoliubov formalism (see Refs. 0, |2jj), which deals 
with coupled equations for the atomic and the molecu- 
lar mean fields, as well as the normal and the anomalous 
densities describing the second-order correlations of the 
non-condensate atomic fields. These correlations are also 
taken into account in the microscopic quantum dynam- 
ics approach used in Ref. Some of these methods, 
however, have no room for incorporating the deactivating 
collisions. The parametric approximation used in Refs. 
[SHI El incorporates both the non-mean-field ef- 
fects and the damping due to deactivating collisions. 

The present work is organized in the following way. 
Section [i] describes the various processes in a hy- 
brid atom-molecule condensate and theoretical methods 
needed for their analysis. The necessary parameters of 
the Cs BEC are estimated in Sec. |n]by a fit of calcula- 
tion results to the experimental data. The formation of 
molecules in the switching scheme, one of the two sweep- 
ing methods used in the experiment Q and discussed in 
Sec. UTI below, is analyzed in Sec. IIIII Optimal conditions 
for the molecular formation are determined in Sec. I1VI 



I. THEORETICAL METHODS 

The effect of Feshbach resonance appears in a BEC of 
Cs atoms when the collision energy of a pair of atoms in 
an open channel is close to the energy of a bound state 
Cs2 (m) in a closed channel (see Ref. jg). The temporary 
formation and dissociation of the resonant (Feshbach) 
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molecular state CS2 (to) can be described as a reversible 
reaction 



Cs + Cs 



CS2 (m) . 



(1) 



The CS2 (to-) is an excited rovibrational state with orbital 
angular momentum I = 4, belonging to an excited state 
of the fine and hyperffne structures. This state can be 
deactivated by an exoergic collision with a third atom of 
the condensate HITollllllTl. 



Cs 2 (to) + Cs^ Cs 2 (d) + Cs 



(2) 



bringing the molecule down to a lower state Cs2 (d) , and 
releasing kinetic energy to the relative motion of the re- 
action products. Although the collision occurs with a 
vanishingly small kinetic energy, rates of such inelastic 
processes remain finite at near-zero energies 0, l25j . A 
variant of this process, involving deactivation by a colli- 
sion with another molecule (rather than an atom), of the 
type 



Cs 2 (to) + Cs 2 (m) -> Cs 2 (d) + Cs 2 (u) 



(3) 



would require a significant molecular density to be effec- 
tive. The two molecular states Cs 2 (d) and Cs 2 (u) can 
be distinct. 

A complete analysis of the processes in a hybrid atom- 
molecule BEC must take into account both the relaxation 
processes due to deactivating collisions and and 
the quantum fluctuations due to dissociation of the res- 
onant molecules to non-condensate atoms QJ. A start- 
ing point of such an analysis can be the quantum equa- 
tion of motion for the atomic field annihilation operator 
ty a (p, t) in the momentum representation. An adiabatic 
elimination of the "dump" states Cs 2 (d) and Cs 2 (u) (see 
Refs. 0, HU) reduces this equation to the Heisenberg- 
Langevin stochastic equation 



ih& a (p, t) = 



2m 2 



+2g*<p m (t) *t (_ P; t ) + iF (p, t) . (4) 



Here m is the atomic mass, e a (t) = — hn(B(t)—Bo) 
is the time-dependent Zeeman shift of the atom in an 
external magnetic field B (t) relative to half the energy 
of the molecular state (which is fixed as the zero-energy 
point), /i is the difference between the magnetic momenta 
of an atomic pair and a molecule, and Bq is the reso- 
nance value of B. The atom-molecule hyperfine coupling 
g is related to the phenomenological resonance strength 
A through \g\ 2 = 27rfi 2 |a a |/xA/TO (see Ref. 0), where 
a a is the background elastic scattering length for atom- 
atom collisions. The deactivation J2J) by atom-molecule 
collisions is represented in Eq. Q by the imaginary term, 
proportional to the deactivation rate k a , as well as by the 
quantum noise source F(p,t), related by a fluctuation- 
dissipation theorem. The quantum noise is required in 
order to maintain the correct commutation relations of 
the atomic field operators. 



In the parametric approximation |l7l |24| . the quan- 
tum fluctuations of the molecular field arc neglected, and 
the molecules are described by a mean field (p m (t). The 
atomic field operator is expressed in this method as 

* a (p, t)=C (t) [i (p, t) yj c (p, t) + it (-p, t) yj s (p, t) 

where the damping factor C (t) takes into account the 
imaginary term in Eq. l|4"|). 



C (t) = exp 



(6) 



and the c-number functions ip CtS (p, t) satisfy the ordinary 
differential equations (with p as a parameter) 



iHipcs (p, t) 



i> c , s {p,t)+2g*(p m (t) # jC (p,t). 

(7) 



The initial conditions ip c {p, to) = 1, ips (P) to) = are in- 
troduced at t = to, assuming the atomic field is then a co- 
herent state of zero kinetic energy. The operators A (p, t) 
can be expressed in terms of the functions ip c ,s (p, t) and 
the quantum noise F (p, t). As a result of ensuing analy- 
sis (see Refs. PJll23); the atomic density comprises the 
sum 



n a (t) = no (t) + n s (t) 
of the densities of condensate atoms 
no (t) = \ip (t)| 2 , 
and of non-condensate (entangled) atoms 



n s (t) = (2irh) 3 / d J pn s (p, t) 



(8) 
(9) 

(10) 



Here 



ip (t) = C (t) [V>c (0, t) ifo (t ) + i> s (0, t) ifl (to)] (11) 

is the atomic condensate mean field. The momentum 
spectrum of the non-condensate atoms 

n s (p, t) = \ijj s (p, t) | 2 [1 + n s (p, t)} + (p, t) \ 2 n s (p, t) 

-2Re ir s (p,t)iPc(p,t) r)c (p,t)U2) 

as well as their anomalous density [encountered in Eq. 
(JTSJl below] 

m s (p, t) = ip s (p, t) ip c (p, t) [1 + 2n s (p, t)} 

-i>l (p, t) n c (p, t) - yj 2 s (p, t) n* c (p, t) , (13) 

are expressed in terms of the auxiliary functions 

t 

n s (p,t) = k a C 2 (t) 



dt' 



to 



r , ( -|^ m (0^(p,OI 2 



(14) 



dt 



Vc (p,t) = k a C 2 (t) J Q^frrn Vma(pJ)P c (P J), 
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which describe the contribution of quantum noise. 

The equation of motion for the molecular mean field 
has the form (see Refs. HEH) 

ihip m (t) = gtfl (t) - i \ -^ n a (*) + k m \tp m (t) | 2 j ip m (t) 



1 



dp [p 2 gm s (p, t) + 2hm\g\ 2 ip m (t)](25) 



where k m is the rate coefficient of the molecule-molecule 
deactivating collisions J2J. The second term under the 
integral over p appears as a result of a renormalization 
procedure (see Refs. 0,0])' necessary in order to regu- 
larize the integral. A numerical solution of Eqs. Q on a 
grid of values of p, combined with Eq. (|15|l . is consistently 
sufficient for elucidating the dynamics of the system. 

The parametric approximation considered above is 
particularly suitable for the analysis of homogeneous sys- 
tems. It can be applied also to inhomogcncous systems 
using a local density approximation, but its application 
to a strongly inhomogeneous expanding BEC meets seri- 
ous difficulties. Fortunately, under proper conditions this 
case can be treated sufficiently well by a mean field ap- 
proach (see Ref. neglecting the atomic field quan- 
tum fluctuations. The applicability of this simpler ap- 
proach can be verified by a comparison of results of 
the parametric and mean-field calculations for the cor- 
responding homogeneous system. 

The expansion of a pure atomic BEC has been con- 
sidered in Ref. j2|| by the introduction of scaled normal 
coordinates 

Pj = r 3 /b, (t) , 1 < j < 3, (16) 

where the scales bj obey the equations 

b j (t) = ^/[b 1 (t)b 2 (t)b 3 (t)b j (t)}, (17) 

in which the uij are the angular frequencies of the har- 
monic trap containing the condensate before expansion. 



The initial conditions bj (t exp ) — 1, bj (t e 



are 



stated at the start of the expansion t cxp . Solutions of 
Eq. 117|) (see Ref. |26j) demonstrate that the expansion 
is ballistic after an acceleration period of ~ min (to J 1 ). 
As shown in Ref. |2^| , the molecules inherit the velocity of 
the atoms they are formed from. The atomic and molec- 
ular mean fields can be represented in terms of rescalcd 
fields $o (Pit) and $ m (p, t), respectively, as 



(fo (r,t)=A(t) $ (p,t)e 



,s 



(18) 



Vm (r,t) = A(t) $ m (p,t) e 2lS , 



where the scaling factor A (t) = (b\ (t) 62 (t) b% (t)) ' 
describes the density reduction, and the phase factor with 



s W = -h2^ r hbj(t) 

.7=1 j 



I J dt'A 2 {t') (19) 



contains most of the contribution of the kinetic energy. 
Here eo = 4^irh 2 a a no/m is a chemical potential of the 
atomic BEC and no is its peak density while the trap is 
on. As a result ( see Ref. (23|), the rescaled mean fields 
obey the set of ordinary differential equations 



ih<S> (p, t) 



e a (t)- -A 2 (t)fc |* m (p,t)| 2 



$0 (p,t) 



-2A(t)g*$*(p,t) $ m (p,t) 



(20) 



ih$ m { Pl t) = -iA 2 (t) 



$ m ( P ,t) + A(t) g$l{p,t). 



The coordinate dependence arises from the use of 
inhomogeneous Thomas-Fermi initial conditions for 
$0 (p, tcxp), while $ m (p, t cxp ) = 0. 



II. PARAMETER ESTIMATES 

Molecules have been formed in the experiments by 
using a very weak Feshbach resonance in Cs situated near 
20 G. The resonance strength and the rate coefficients 
of atom-molecule and molecule-molecule deactivation are 
unknown and are estimated here by a fit of the calculation 
results to the experimental data. 

The magnetic field has been varied in these experi- 
ments in two manners. In the ramping scheme the mag- 
netic field has been swept through resonance with a fixed 
ramp speed. In the switching scheme the magnetic field 
has been tuned to a value -Btest in the vicinity of the res- 
onance and then held for a fixed time fholdj starting from 
the value B sta , T t = B ±0.5 G for Btest ^ -Bo, respectively, 
so that the resonance should not be crossed. Due to finite 
response time the magnetic field variation is represented 
by an exponential function, 



B (i) = B t , 



{B s 



- Btest) exp [(t - t cxp ) /1.54 ms] . 

The switching scheme has been applied both to the 
trapped and the expanding BEC. In the last case, the 
magnetic field variation has been started from -B sta rt si- 
multaneously with the expansion at t = t exp . 

Consider first the switching scheme for a trapped BEC, 
resulting in a condensate loss with a negligible molecular 
formation. This experiment is similar to the slow-sweep 
Na experiments [27| and the 87 Rb experiments [2^], in 
which the field was stopped short of resonance, too. In 
those cases the condensate loss is determined by simple 
analytical expressions involving the product of the res- 
onance strength A and the atom-molecule deactivation 
rate coefficient k a (see Refs. 0,0,|2{J). In the present 
case, these analytical expressions are inapplicable owing 
to the non-linear magnetic field variation. However, nu- 
merical calculations still demonstrate a dependence on 
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FIG. 1: Switching-scheme atom loss fraction calculated for 
the trapped BEC with A = 8 mG, k m = 1.5 x 10" 9 cm 3 /s, 
k a = 1 x 10" 11 cm 3 /s (solid line) and k a = 2 x 10~ n cm 3 /s 
(dot-dashed line). The results for k m = 1 x 10 -10 cm 3 /s are 
calculated with A = 8 mG, k a — 1 x 10 -11 cm 3 /s (dashed 
line) and A = 2 mG, k a = 4 x 1CT 11 cm 3 /s (pluses). The 
open and solid circles represent the experimental data of Ref. 
U measured below and above the resonance, respectively. 



FIG. 2: Ramping-scheme atom loss fraction (1) and conversion 
efficiency (c) calculated for the expanding BEC with k m = 
1.5 x 10~ 9 cm 3 /s, k a A = 8 x 10 -11 mG cm 3 /s, and A = 8 
mG (solid lines), A = 7 mG (dashed lines), or A = 9 mG 
(dot-dashed lines). The results for A = 8 mG and k a = 
1 x 10 -11 cm 3 /s are represented by pluses (k m — 2 x 10~ 9 
cm 3 /s) and crosses (fc m = 1 x 10 -9 cm 3 /s). The open and 
solid circles represent the experimental data of Ref. jjj for loss 
and conversion, respectively. 



the product k a A provided moderate molecule- molecule 
deactivation rates are assumed. (This is demonstrated 
by the dashed line and pluses in Fig.^). The comparison 
with the experimental data leads to an estimated value 
of k a A w 8 x 10 -11 mG cm 3 /s. Much higher values of 
k m do not lead to such a good fit. Although the conden- 
sate loss becomes dependent on a variation of k a and A, 
keeping the product k a A fixed, this variation does not 
improve the fit. 

Consider now the ramping scheme. It has been applied 
in experiments to the expanding BEC, measuring both 
the numbers of remained atoms and formed molecules. 
In the fast-decay approximation (see Refs. [lUlS)) the 
atomic condensate loss due to the resonance crossing is 
determined by the resonance strength A only and is in- 
dependent of the deactivation rates. Although the an- 
alytical expressions of Refs. 10, :12j are inapplicable to 
the case of an expanding BEC, the results of numerical 
calculations demonstrate a low sensitivity of the BEC 
loss to the deactivation rates. A comparison with the 
experimental data leads to the estimated value of A « 8 
mG (see the upper graphs in Fig. EJ) . Together with the 
above estimate for k a A this leads to k a w 1 x 10~ n 
cm 3 /s. The experimental data points presented in Fig. [21 
were obtained with a magnetic field ramp that has been 
started at t cxp , using a magnetic field value such that 
the resonance is crossed in 10 ms and the populations 
are measured 20 ms after t exp [30|. 

The remaining unknown parameter, the molecule- 
molecule deactivation rate coefficient k m , can be esti- 
mated by a comparison of the calculation results with 



the experimental data for the number of atoms converted 
to molecules (see lower graphs in Fig. EJ). This leads to 
the estimated value of k m w 1.5 x 10~ 9 cm 3 /s. The 
graphs in Fig. |2 are practically insensitive to k a because 
of its much lower value. A look back at Fig. ^ f° r the 
switching scheme shows that, at most, k a is bounded by 
k a < 2 x 10 -11 cm 3 /s, when the far-off-resonance results 
are regarded. In this range, the results are not much sen- 
sitive to the value of k m . However, nearer resonance, the 
higher value of k m would lead to an overestimate of the 
loss, independently of k a and A. 

Let us compare these values with what is known re- 
garding other atoms. The estimated value of the rate 
coefficient of atom-molecule deactivation k a w 10~ u 
cm 3 /s is several times lower than the corresponding val- 
ues of 5.5 x 10 -11 cm 3 /s for Na measured in Ref. [g 
and 7 x 10 -11 cm 3 /s for 87 Rb estimated in Ref. |29j . 
A molecule-molecule deactivation rate coefficient k m ~ 
1.5 x 10~ 9 cm 3 /s exceeds the corresponding value for Na 
2.5 x 10~ n cm 3 /s measured in Ref. 6| by two orders of 
magnitude. These differences may be related to the large 
orbital angular momentum (I — 4) of the resonant molec- 
ular state in the present Cs case, while for Na and 87 Rb 
/ = 0. But there may be another reason. 

All the estimates made here concerning the data for 
molecular conversion are based on the suggestion that the 
resonant beam, blasting out the atoms in the experiments 
fH, does not affect the molecules. (This procedure was 
used to separate the atoms from molecules.) If some part 
of the molecular population is removed by the blasting 
pulse, the experimental results can be explained using 
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FIG. 3: Time dependence of the atomic (a) and molecular 
(m) condensate populations (scaled to the initial atomic one) 
calculated in a ramping scheme for a homogeneous BEC with 
the initial atomic density 6.4 x 10 11 cm -3 , ramp speed 3 G/s 
(in the backward sweep), A = 8 mG, k a = 1 x 10~ n cm 3 /s, 
and k m — 1 x 10~ 9 cm 3 /s with the parametric (solid lines) and 
mean-field (crosses) approximations. The dashed lines repre- 
sent the mean field results for the expanding BEC with the 
initial atomic density 7.6 x 10 13 cm -3 . The non-condensate 
atom population calculated with the parametric approxima- 
tion is plotted by the dot-dashed line. 



lower values of k m . 

All the calculations above were performed with the 
mean field approximation. Figure [3] compares results of 
the parametric and mean-field calculations for the initial 
atomic density of 6.4 x 10 11 cm' 3 , corresponding to the 
mean density at the resonance crossing for the slowest 
ramp speed of 3 G/s used in the experiments 2\. This 
figure demonstrates that the temporary non-condensate 
atom population persists only about 1 ms and after this 
short time the results of parametric and mean-field calcu- 
lations for the homogeneous case coincide within a good 
accuracy. The expansion reduces the deactivation losses 
compared to the ones for the homogeneous case. 



III. MOLECULAR FORMATION IN THE 
SWITCHING SCHEME 

The highest conversion efficiency, reaching beyond 
30%, has been observed in the experiments for the 
switching scheme in an expanding BEC. The character- 
istic time of atom-molecule relaxation t am is determined 
by the coupling terms in Eq. (|2"0)l as 



\Ag$ \ 



27r|a a |^A 



-no 



1/2 



(22) 



where n (t) = \ip (t) | 2 is non-rescaled atomic conden- 
sate density. Even for hq — 10 12 cm~ 3 , corresponding 



to the resonance approach time in the experiments 0], 
tam ~ 1 m s is less than the characteristic time of the 
magnetic field variation (1.54 ms). Therefore the evolu- 
tion of the atom-molecule condensate is adiabatic. Ne- 
glecting the deactivation, it can be described by a quasi- 
stationary solution of the coupled Gross-Pitaevskii equa- 
tions (see Ref. 8j) 



n 



36 



24 - e z + 24 



( 2e 2 ~ 2ey / e 2 + 24 + 24 
144 V 



(23) 



corresponding to the pure atomic BEC (n m = 0) above 
the resonance at e — ► oo. Here n m = |vm| 2 is the non- 
rescaled molecular density, n = no + 2n m is the total 
density of atoms, and 



2e„ 



\g\Vn 



Bo — B t 
W 



SB 



2-irh 2 



-no 



1/2 



(24) 

is a dimensionless detuning (generally time-dependent). 
The later fast sweep of the magnetic field going under 
the resonance, used in the switching scheme, conserves 
the atomic and molecular densities (|23|l acquired at B = 

■Btcst ■ 

The conversion efficiency reaches its maximum 



(25) 



when the magnetic field levels off at e = 0. In the case of 
large detunings \B tcs t — B \ ^> SB the molecular density 
decreases as n m ~ nSB 2 / (B tC st — ^o) 2 . Therefore the 
substantial molecular population and condensate losses 
induced by the deactivating collisions can take place only 
while latest - B \ < SB w 0.2 mG for n = 10 12 cm" 3 . 

This conclusion is confirmed by the results of the nu- 
merical calculations (see Fig. 0}. These results, how- 
ever, predict a maximal conversion efficiency of 36% for 
B tcst = Bo — 0.15 mG, when the resonance is crossed. 
The smooth time dependence of the atomic and molecu- 
lar populations (see Fig.0 is in agreement with the adi- 
abatic evolution mentioned above. On decreasing further 
Btcst the conversion efficiency decreases, demonstrating 
Rabbi oscillations due to non-adiabatic effects. The re- 
sults in Fig. 0] correspond to a resonance approach from 
above {B start = B + 0.5 G). In the mean field approxi- 
mation used here, neglecting molecular dissociation with 
formation of non-condensate atoms, the same results are 
reached by approaching the resonance from below. 

Although the theory describes the peak conversion effi- 
ciency and condensate losses observed in the experiments 
|2j, the actual experimental width of the resonance in 
conversion and in loss, of 2 mG, is about an order of mag- 
nitude more than the theoretical one. This disagreement 
can be related to magnetic field variation mentioned in 
Ref. 0. 
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FIG. 4: (a) Atom loss fraction (dashed line) and conversion 
efficiency (solid line) calculated as functions of the resonance 
detuning for the switching scheme in the expanding BEC with 
k m = 1.5 x 10~ 9 cm 3 /s, k a = 1 x 1CT 11 cm 3 /s, and A = 8 
mG. The open and solid circles represent the experimental 
data of Ref. for loss and conversion, respectively. The re- 
sults of calculations for a rippled magnetic field with the phase 
X = 1-25tt are plotted by dot-dashed (conversion) and dashed 
(loss) lines, (b) The loss (dashed line) and conversion (solid 
line) for x — 1.85ir. The results for inhomogeneous magnetic 
field with gradient 1 G/cm are plotted by dot-dashed (con- 
version) and dashed (loss) lines. In both parts the triangles 
and circles represent the experimental data of Ref. |2j for loss 
and conversion, respectively. 











0.8 " 






s o c 
/ o ° o 


c0.6 " 
o 

o 
cd 

~0.4 " 


/ 

/ 

/ 

/ 

/ 

/ 


y 


A , 


0.2 


/ 

/ / 
/ / 
/ / 
// 
// 

. i. » * . 


o 
• 


■i 

:i 
■i 

.7 



10 15 

thold < ms ) 



FIG. 5: Time dependence of the atom loss fraction (dashed 
line) and conversion efficiency (solid line) calculated for the 
switching scheme in the expanding BEC with A = 8 mG, 
k m = 1.5 x 10" 9 cm 3 /s, k a = 1 x 10~ n cm 3 /s, and S tes t = 
—0.15 mG. The open and solid circles represent the experi- 
mental data of Ref. 2] for loss and conversion, respectively. 
The results of calculations for a rippled magnetic field with 
the phase x — 1-25-7T and B tas t = —0.4 mG are plotted by 
dot-dashed (conversion) and dashed (loss) lines. 



FIG. 6: Time dependence of the rippled magnetic field with 
the phases x — 1-257T, x = 1-85tt and without ripples for 
Bstart = 0.5 G (solid lines) and Bstart = —0.5 G (dashed 
lines). An arbitrary value of Bo is shown here for reference 
by the dotted line. 

An ambient magnetic field ripple with a frequency of 50 
Hz and an amplitude of 4 mG leads to a quite perceptible 
effect. Although the experiment has been synchronized 
with the ripple, the results of calculations demonstrate a 
strong dependence on the ripple phase. In these calcula- 
tions the magnetic field time dependence has the form 

B (t) = Btcst + (-B s tart - Btest) exp [(t - t cxp ) /1.54 ms] 

+2 mG sin [2?r x 50 Hz(t - i cxp ) + x](26) 

The variation of the phase \ leads to a shift of the peak 
and changes the shape of the magnetic field dependence 
of the loss and molecule formation (see Fig.0J). 

The ripple may lead to the formation of molecules due 
to the resonance crossing in a backward direction (see 
Fig-©- The crossing can be non-adiabatic, leading to an 
oscillating time dependence of the atomic and molecular 
populations (see Figs.|3]and[7I). The oscillations are sensi- 
tive to the ripple phase x as well. For some values of x, e. 
g. x = 1-85, the resonance can be crossed a second time, 
but in the forward direction, even though the first cross- 
ing occurred by approaching from above (B stalt > Bq), 
given the arbitrary value of Bq shown in Fig. |SJ In this 
case, a rather long hold-on-time between the crossing and 
measurement leads to sharp Rabbi oscillations. These os- 
cillations, however, can be averaged by a magnetic field 
inhomogeneity with a rather small gradient of l mG/cm 
(see Figs. 0] and EJ). Effects of both oscillations and in- 
homogeneity can broaden the resonance to about l mG, 
which is still less than the experimentally observed value 
of 2 mG. The additional broadening can be related to 
the uncontrolled magnetic field variations of about l mG 
mentioned in Ref. [2], the behavior of which is unclear. 

The magnetic-field time dependence Il2(it can lead to 
different results in approaching the resonance from below, 
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hold 1 



FIG. 7: Time dependence of the atom loss fraction (dashed 
line) and conversion efficiency (solid line) calculated for a 
rippled magnetic field with the phase x = 1.857T, and with 
St«t = 1.7 mG, k m = 1.5 xl(T 9 cm 3 /s, k a = lxlO" 11 cm 3 /s, 
and A = 8 mG. The results for an inhomogeneous magnetic 
field with gradient 1 G/cm are plotted by dot-dashed (con- 
version) and dashed (loss) lines. 



with Bstart < Bq i due to an additional resonance crossing 
occurring earlier. This forward crossing can lead to ad- 
ditional condensate losses due to molecular dissociation 
into pairs of non-condensate atoms. This dissociation 
can, however, be reversible, as the non-condensate atoms 
can associate to molecules during the following backward 
crossing. A correct analysis of these effects requires the 
application of non-mean-field calculations to the expand- 
ing BEC, in order to obtain reliable conversion estimates. 



IV. OPTIMAL CONDITIONS 

Having a model consistent with the experimental data, 
we can proceed to determine optimal conditions for the 
molecular formation. In the switching scheme, when 
the magnetic field adiabatically approaches the resonance 
and further suddenly crosses it, the conversion efficiency 
is restricted by the value of ^ (see Eq. (|25H . An adiabatic 
crossing, as in the ramping scheme, allows for higher re- 
sults since the adiabatic state (|23|l corresponds to a total 
conversion (n m = n/2) in the limit of e — > oo below the 
resonance. As in the cases of Na [13] and 87 Rb p5l2^ . 
the optimal ramp speed is determined by a balance of 
the atomic association, decreasing at faster sweeps, and 
deactivation losses, increasing at slower sweeps. Thus, 
a sudden sweep would lead to no association, while an 
infinitely slow adiabatic sweep would lead to total asso- 
ciation, accompanied by a total deactivation loss during 
the infinite time. The optimal density is determined by 
a concurrence of deactivation losses, increasing at high 
densities, and dissociation into non-condensate atoms, in- 




initial atomic density (cm 



FIG. 8: (a) Conversion efficiency (solid line) at the optimal 
ramp speed (dashed line), (b) The lifetime of the molecular 
condensate (solid line) and the time after the resonance cross- 
ing when the peak molecular density is reached (dashed line) , 
using the appropriate optimal ramp speed. 
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FIG. 9: Conversion efficiency (solid line) and the lifetime of 
the molecular condensate (dashed line) as a function of the 
ramp speed for the initial atomic density 5 x 10 10 cm -3 . 



creasing at low densities. An analysis taking account of 
the latter loss process requires a non-mean-field theory, 
such as the parametric approximation. It was mentioned 
in Sec.Q]above that this approximation has its limitations 
in dealing with expanding gases. However, the effect of 
increasing density due to expansion is less important at 
the rather low optimal density, stated below, compared 
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time (ms) 

FIG. 10: Time dependence of the atomic (dashed line) and 
molecular (solid line) condensate populations (scaled to the 
initial atomic one) calculated for a homogeneous BEC with 
the initial atomic density 5 x 10 10 cm" 3 and ramp speed 35 
mG/s (in the backward sweep). The non-condensate atom 
population is plotted by the dot-dashed line. The resonance 
is crossed at t = 0. 

to the conditions pertaining to Fig. [3] For these reasons 
the optimal conditions are determined by using the para- 
metric approximation for a homogeneous non-expanding 
BEC. 

The results presented in Fig. |H1 show an optimal initial 
atomic density of 5 x 10 10 cm -3 . This value corresponds 
to the mean density of a Thomas- Fermi distribution with 
a peak density of 1.25 x 10 11 cm" 3 . Under the conditions 
of the experiments 0, this density can be reached after 
35 ms of expansion. The optimal ramp speed of 35 mG /s 



is much slower than the one used in the experiments. The 
high conversion efficiency does not change much when 
the ramp speed is varied by about a order of magnitude 
(see Fig.[|JJ). Figure ITTJ1 shows that the molecular density 
reaches its maximum in about 8 ms after the resonance 
crossing, or about 0.25 mG below the resonance, and 
the ramp should be started about 2 mG above it. This 
figure demonstrates also a substantial population of non- 
condensate atoms, justifying the necessity to use non- 
mcan-ficld calculations. 



Conclusions 

The loss of Cs atomic BEC and formation of a molecu- 
lar BEC observed in the experiments can be described 
by a mean-field theory of an expanding atom-molecule 
BEC. A fit of the calculation results to the experimental 
data leads to estimated values of the resonance strength 
and rate coefficients for atom-molecule and molecule- 
molecule deactivating collisions. At small detunings the 
results are sensitive to magnetic field ripple and inhomo- 
geneity. A determination of optimal conditions for the 
molecular formation requires a non-mean-field paramet- 
ric approximation, taking into account the dissociation 
of molecules into non-condensate atomic pairs. A con- 
version efficiency of 60% is predicted for lower densities 
and slower sweeps than the ones used in the experiments. 



Acknowledgments 

The authors are very grateful to R. Grimm, H.-C. 
Nagerl, and C. Chin for helpful discussions and clarifi- 
cation of experimental details. 



[1] J. Herbig, T. Kraemer, M. Mark, T. Weber, C. Chin, 

H.-C. Nagerl, and R. Grimm, Science 301, 1510 (2003). 
[2] M. Mark, T. Kraemer, J. Herbig, C. Chin, H.-C. Nagerl, 

and R. Grimm, Europhys. Lett. 69, 706 (2005). 
[3] S. Diirr, T. Volz, A. Marte, and G. Rempe, Phys. Rev. 

Lett. 92, 020406 (2004). 
[4] S. Diirr, T. Volz, and G. Rempe, Phys. Rev. A 70, 

03160RR) (2004). 
[5] K. Xu, T. Mukaiyama, J. R. Abo-Shaeer, J. K. Chin, D. 

E. Miller, and W. Ketterle, Phys. Rev. Lett. 91, 210402 

(2003). 

[6] T. Mukaiyama, J. R. Abo-Shaeer, K. Xu, J. K. Chin, and 
W. Ketterle, Phys. Rev. Lett. 92, 180402 (2004). 

[7] C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Na- 
ture 424, 47 (2003); M. Greiner, CA. Regal, and D. S. 
Jin, Nature 426, 537 (2003); K. E. Strecker, G. B. Par- 
tridge, and R. G. Hulet, Phys. Rev. Lett. 91, 080406 
(2003); J. Cubizolles, T. Bourdel, S. J. J. M. F. Kokkel- 
mans, G. V. Shlyapnikov, and C. Salomon, Phys. Rev. 
Lett. 91, 240401 (2003); T. Bourdel, L. Khaykovich, J. 
Cubizolles, J. Zhang, F. Chevy, M. Teichmann, L. Tar- 



ruell, S. J. J. M. F. Kokkelmans, and C. Salomon,, Phys. 
Rev. Lett. 93, 050401 (2004); S. Jochim, M. Barten- 
stein, A. Altmeyer, G. Hendl, C. Chin, J. H. Denschlag, 
and R. Grimm, Phys. Rev. Lett. 91,240402 (2003); S. 
Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, S. Riedl, 
C. Chin, J. H. Denschlag, and R. Grimm, Science 302, 
2101 (2003); M.W. Zwierlein, C. A. Stan, C. H. Schunck, 
S.M. F. Raupach, S. Gupta, Z. Hadzibabic, and W. Ket- 
terle, Phys. Rev. Lett. 91, 250401 (2003); J. Kinast, S. 
L. Hemmer, M. E. Gehm, A. Turlapov, J. E. Thomas, 
Phys. Rev. Lett. 92, 150402 (2004). 
[8] E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- 

man, Phys. Rep. 315, 199 (1999). 
[9] F. H. Mies, E. Tiesinga, and P. S. Julienne, Phys. Rev. 
A 61, 022721 (2000). 

[10] V. A. Yurovsky, A. Ben-Reuven, P. S. Julienne and C. J. 
Williams, Phys. Rev. A 60, R765 (1999); 

[11] F. A. van Abeelen, and B. J. Verhaar, Phys. Rev. Lett. 
83, 1550 (1999). 

[12] V. A. Yurovsky, A. Ben-Reuven, P. S. Julienne and C. J. 
Williams, Phys. Rev. A 62, 043605 (2000). 



9 



[13] M. Shapiro and P. Brumer Principles of the Quantum 

Control of Molecular Processes (John Wiley and Sons, 

New York, 2003). 
[14] C. P. Koch, J. P. Palao, R. Kosloff, and F. Masnou- 

Seeuws, Phys. Rev. A 70, 013402 (2004); E. Luc-Koenig, 

R. Kosloff, F. Masnou-Seeuws, and M. Vatasescu, Phys. 

Rev. A 70, 033414 (2004). 
[15] D. S. Petrov, Phys. Rev. A 67, 010703(R) (2003); D. S. 

Petrov, C. Salomon, and G.V. Shlyapnikov, Phys. Rev. 

Lett. 93, 090404 (2004); Phys. Rev. A 71, 012708 (2005). 
[16] E. Pazy, A. Vardi, Y. B. Band, Phys. Rev. Lett. 93, 

120409 (2004); J. Chwedehczuk, K. Goral, T. Kohler, 

P. S. Julienne, Phys. Rev. Lett. 93, 2 60403 (2004); I. 

Tikhonenkov, A. Vardi, cond-mat/0407424 
[17] V. A. Yurovsky and A. Ben-Reuven, Phys. Rev. A 67, 

043611 (2003). 

[18] U. V. Poulsen and K. Molmer, Phys. Rev 63, 023604 

(2001) ; J. J. Hope, M. K. Olsen, and L. I. Plimak, Phys. 
Rev. A, 63, 043603 (2001); J. J. Hope, and M. K. Olsen, 
Phys. Rev. Lett., 86, 3220 (2001); M. K. Olsen, L. I. Pli- 
mak, and M. J. Collett, Phys. Rev. A 64, 063601 (2001); 
K. V. Kheruntsyan and P. D. Drummond, Phys. Rev. 
A, 66, 031602(R) (2002); M. K. Olsen and L. I. Plimak, 
Phys. Rev. A 68,031603(R) (2003); M. K. Olsen, Phys. 
Rev. A 69, 013601 (2004); M. K. Olsen, A. S. Bradley, 
and S. B. Cavalcanti, Phys. Rev. A 70, 033611 (2004); 
K. V. Kheruntsyan, M. K. Olsen, and P. D. Drummond, 
|cond-mat /0407363 

[19] M. Holland, J. Park, and R. Walser, Phys. Rev. Lett. 86 
1915 (2001). 

[20] S. J. J. M. F. Kokkelmans, M. J. Holland, Phys. Rev. 

Lett. 89 180401 (2002). 
[21] T. Kohler and K. Burnett, Phys. Rev. A 65, 033601 

(2002) ; T. Kohler, T. Gasenzer and K. Burnett, Phys. 



Rev. A 67, 013601 (2003); K. Goral, T. Kohler, S. A. 
Gardiner, E. Tiesinga, and P. S. Julienne, J. Phys. B 37, 
3457 (2004); K. Goral, T. Kohler, and K. Burnett, Phys. 
Rev. A 71, 023603 (2005). 
[22] A. Vardi, V. A. Yurovsky, and J. R. Anglin, Phys. Rev. 
A 64, 063611 (2001); V. A. Yurovsky, Phys. Rev. A 
65, 033605 (2002); V. A. Yurovsky, A. Ben-Reuven, 
and P. S. Julienne, Phys. Rev. A 65, 043607 (2002); V. 
A. Yurovsky and A. Ben-Reuven, J. Phys. B 36, L335 
(2003). 

[23] V. A. Yurovsky and A. Ben-Reuven, Phys. Rev. A 70, 
013613. (2004). 

[24] V. A. Yurovsky, in Soft Condensed Matter: New Research 
(Nova Science Publishers, Hauppauge, in press). 

[25] R. C. Forrey, V. Kharchenko, N. Balakrishnan and A. 
Dalgarno, Phys. Rev. A 59, 2146 (1999); P. Soldan, M. 
T. Cvitas, J. M. Hutson, P. Honvault, and J.-M. Launay, 
Phys. Rev. Lett. 89, 153201 (2002). 

[26] Yu. Kagan, E. L. Surkov and G. V. Shlyapnikov, Phys. 
Rev. A 54, R1753 (1996); Y. Castin and R. Dum, Phys. 
Rev. Lett. 77, 5315 (1996). 

[27] S. Inouye, M.R. Andrews, J. Stenger, H.-J. Miesner, D. 
M. Stamper-Kurn, and W. Ketterle, Nature 392, 151 
(1998); J. Stenger, S. Inouye, M.R. Andrews, H.-J. Mies- 
ner, D.M. Stamper-Kurn, and W. Ketterle, Phys. Rev. 
Lett. 82, 2422 (1999). 

[28] A. Marte, T. Volz, J. Schuster, S. Diirr, G. Rempe, E. 
G. M. van Kempen, and B. J. Verhaar, Phys. Rev. Lett. 
89 283202 (2002). 

[29] V. A. Yurovsky and A. Ben-Reuven, Phys. Rev. A, 67, 
050701 (R) ( 2003). 

[30] C. Chin, private communication. 



